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Abstract 

We propose a precise definition of multidimensional fluids generated by self-gravitating ex- 
tended objects such as strings and membranes: a p-dimensional perfect fluid is a smooth involutive 
p-dimensional distribution on a spacetime, each integral manifold of which is a timelike, connected, 
immersed submanifold of dimension, p - representing the history of a {p — l)-dimensional ex- 
tended object. This geometric formulation of perfect fluids of higher dimensions naturally leads to 
the associated stress-energy tensor. Furthermore, the laws of temporal evolution and symmetries 
of such systems are derived, in general, from the Einstein field equations and the integrability 
conditions. We also present a matter model based on a 2-dimensional involutive distribution, 
and it is shown that the stress-energy tensor for self-gravitating strings gives rise to a non-trivial 
spherically symmetric spacetime with a naked singularity. 



1 Introduction 



The purpose of this work is to develop a general relativistic theory of multidimensional 
fluids as sources of spacetime curvature. The basic ingredients of such a fluid are (p — 1)- 
dimensional spatially extended objects called p-branes - where p = 1, 2, 3 correspond to point 
particles, strings and membranes respectively. More precisely, a p-brane is a timelike, con- 
nected, p-dimensional C°°-manifold immersed in a spacetime - representing the temporal 
evolution of a (p — 1) -dimensional extended object. Then, for a fixed p, a multidimensional 
fluid is defined by a smooth involutive p-dimensional distribution on a spacetime, each inte- 
gral manifold of which is a p-brane. Thus, a multidimensional fluid naturally generalizes the 
model of a coUisionless gas of point particles by a congruence of world hues (1-dimensional 
distribution). [In this paper, all manifolds, tensor fields on them, and all maps from one 
manifold to another will be C°°. Also, we define a spacetime as a non-compact, connected, 
oriented and time-oriented n-dimensional manifold, M, endowed with a Lorentz metric, g. 
For n — A, the corresponding spacetime will be denoted (M^, g)] 

In a class of field theoretic models^""^ of the early universe, p-branes appear as 'topolog- 
ical defects' with characteristic rest-mass per unit (p — 1) -dimensional spatial volume. The 
existence of such extended objects could be a possible source of density perturbations, and 
hence may provide a causal mechanism for generating the observed large scale structure of 
the universe^. Thus any scheme, based on general relativity, for investigating the role of mul- 
tidimensional fluids in the evolution of the universe, requires specification of a stress tensor, 
solutions of Einstein's field equations, description of the behaviour of other forms of matter 
in the vicinity of the extended objects, and characterization of the resulting spacetimes and 
their singularities. 

In order to carry out this program, we need a precise form of the fluid stress tensor, T, 
which is formally a symmetric (0, 2) -tensor field on a spacetime {M,g). Physically T replaces 
and unifies the concepts of energy density, momentum density, energy flux and momentum 
flux. These quantities are observer dependent. [An ohservei is a future-pointing timelike 
curve 7 : / — > M (7 C R is an open interval) such that Vs e 7, the tangent vector 
7*s e T^(^g)M satisfies g{'y*s,l*s) — ~1- An instantaneous observer {x,Z) at x e M, is a 
future-pointing timelike unit vector Z e T^M.] When an instantaneous observer {x, Z) in 
{M^,g) measures, for instance, the energy density in any unit 3- volume of the local rest 
space Z-^ = {X e TxM'^\g{X, Z) — 0}, T(Z, Z) corresponds to the measured energy density. 
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Also, for all known forms of matter T{Z, Z) > 0, V instantaneous observer {x, Z) (and hence 
by continuity T(X,X) > 0, V causal X e T^M^) Vx e M^. This operational definition 
uniquely specifies T in the following sense^: 

Theorem : If the symmetric (0,2) tensor fields T and T' on a spacetime {M,g) satisfy 
T(Z, Z) — T'{Z, Z) for all instantaneous observers {x, Z) then T = T'. 

Hence measured energy density T{Z,Z) naturally motivates the following 

Definition 1 : A stress-energy tensor on spacetime M is a symmetric (0, 2)-tensor field T on 
M such that J{X, X) > for all causal X e T^M, Va; e M. 

Based on the Theorem above, we shall motivate the definition, (1.3), of the stress tensor 
for a collisionless gas of point particles (of mass m) in a way that is suitable for gener- 
alization to multidimensional fluids. Such a fluid on {M^,g) is a congruence of integral 
curves of a nowhere vanishing (energy- momentum) vector held, P, on a spacetime region - 
where g{P, P) = — m^. This conflguration of (integral) curves is a 1-dimensional involutive 
distribution on (or 1-foliation of M^). 

The above geometric structure of a particle flow suggests that a collection of non-coUiding 
extended particles (non-intersecting p-branes) in a spacetime, (M", g), can be modelled by a 
p-dimensional {p > 1) foliation of {M'",g) - where the timelike integral manifolds represent 
the p-branes. After a brief introduction to foliations in section 2, we then characterise a 
multidimensional fluid, in section 3, as a p-foliation determined locally by a nowhere zero 
decomposable p-form u = o"(Vi A ■ ■ ■ A Vp)/fi where {Va}^=;^ are metric dual of local vector 
flelds {Va} giving (local) bases for tangent spaces of integral manifolds of the j9-foliation, 
det[5f(Va, Vf,)] = — /i^ 7^ 0, and a is the characteristic rest mass per unit {p — l)-dimensional 
spatial volume of the p-branes satisfying Q{uj,uj) = — ^ 0. Here, Q, is a scalar product on 
the vector space of differential p-forms, A^^M"') [see Appendix]. 

This local fohation p-form, together with a density function r] [(3.1)] allows us to 
specify, locally, in a smooth way the number of p-branes of the fluid system in spacelike 

sections of a spacetime. Such a description of a multidimensional fluid in terms of (uj, rj) 
naturally leads to the associated stress tensor [(3.15)] with the following local representation: 
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where, i^, is the interior contraction operator on differential forms with respect to any local 
basis vector fields {Xa} with the corresponding dual basis {e"}. 

In section 4 we consider such a stress tensor as a possible source of spacetime curvature, 
and derive its dynamical consequences from the Einstein field equations. In particular we 
have shown that the foliation p-form u; satisfies a 'conservation law' [Proposition 4] : 

S{ricu) = 

and each integral submanifold (p-brane) determined by uj has vanishing mean curvature, 
H = [Proposition 3]. It is also shown that if the spacetime admits a Killing vector field 
K then the world density, rj as well as uj are invariant with respect to the local isometry 
generated by K [Proposition 1 and Proposition 2] : 

IkT] = ; Lko; = 

[Here S and L are coderivative and Lie-derivative operators respectively]. 

These properties can be used to solve for the foliation p-form, uj, and the world density, 
f], which specify the multidimensional fluid as well as the spacetime metric. In section 5, 
we have considered a matter model based on a 2-foliation where self-gravitating extended 
particles do indeed give rise to a non-trivial spacetime with a naked singularity. A new class 
of gravitational collapse problems is also presented. 

1-Dimensional Perfect Fluids 

In this subsection we motivate the definition of the stress tensor associated with an 1- 
dimensional perfect fluid (flow of point particles) on {M'^,g), from the viewpoint of the 
uniqueness Theorem stated above. For m e [0, oo), a particle of mass m is a future-pointing 
curve 7 : / — such that g{'y*s:7*s) — —rn^-, Vs e /. Here m 7^ is the analogue 
of Newtonian inertial mass and m = is allowed. The vector field, 7*, over 7 is called 
energy-momentum of the particle. Then for an instantaneous observer (7(5), Z), we have 
the orthogonal decomposition of 7*5 G Xy(s)M*^ : 

7*, = eZ + p (1.1) 

where e = —g(j^,,Z) > is the energy and p e Z-^ is the momentum of the particle as 
measured by (x, Z), and hence the Newtonian velocity of a particle with respect to Z is 
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given by V = p/e e Z-^. Now, if we have enormous number of particles, each having the 
same mass m e [0, oo) and the energy-momenta, then we may describe such a system on 
by the following 

Definition 2 : An 1-dimensional perfect fluid (P, 77, m) on [M^^g) consists of a function 77 : 
— > [0, 00) called wovld density and an energy-momentum vector field P : — > TM^ 
such that each integral curve of P is a particle of mass m, and the integral of the number 
density 3-form 

n = ★(77P) (1.2) 

over a spacelike section C defines the total number of particles in V^. Associated 
with (P,?7, m) is the stress-energy tensor 

T^rjP^P (1.3) 

[Here ★ is the Hodge operator induced by the metric on M^, and P is metric dual of P.] 

Remark : It follows from Definition 2 and the definition of a particle (of mass m) that P 
is future-pointing and g{P, P) = —m^. Motivation for (1.3) : T is symmetric, smooth and 
for all X G Tj-M^, Tx{X,X) = rix[g{P, X)]'^ > 0. Thus T is a stress-energy tensor, by 
Definition 1. We now explain in what sense the measured energy density is Tx{Z,Z) for 
every instantaneous observer {x, Z) e T^M^. 

Given an observer {x,Z), (1.2), by regarding part of T^M"^ (using exponential map) as 
a part of for a sufficiently small neighbourhood of a; G M'^, where curvature tensor is 
negligible^. Then, given a set of linearly independent vectors Xi, X2, X3 G Z-^ (rest space 
of Z), the world density function rj can be interpreted as follows : the number of particles 
measured by {x, Z) in the parallelopiped [X1X2X3] C Z-^ C T^M"^ is the number of integral 
cuves of P crossing the parallelopiped [X1X2X3], and is (approximately) given by 

|n(Xi, X2, Xs)! = VxMP, Xi, X2, Xs)! 

= Vxe\n{Z, Xi, X2, Xs)! (1.4) 

where = ^1 is the volume form on M^, \fl{Z, Xi, X2, X^)] is the 3-volume of the 
parallelopiped [X1X2X3], and from the Definition 2 and (1.1), P is given by the following 
orthogonal decomposition 

P, = eZ + p (1.5) 
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An alternative way to calculate the particle number density in is to project, n^; [(1.2)], 
into !s?{Z^) - the vector space of 3-forms on - by the R-hnear map, 
Hz : A=^(7;m4) — > A^(^^), where = 1 + Z A t^. By (1.5), 

^z = ^zV{r]Hx = r]x^^^ (1.6) 

Then the number of particles (x, Z) measures in any unit volume of the local rest space Z^^ 
is given by 

\\iVz\\ = \Q{^z.^z)\^ ^y]x^ (1.7) 

where Q [see Appendix] is the non-degenerate symmetric bilinear form, induced by on the 
vector space Ap(M'^) of differential p-forms. From (1.5) and (1.7) we compute the energy 
density U measured by (x, Z) : 

U = Me (1.8) 
Thus T is a stress tensor for the particle flow (P,77,m), and by (1.3), (1.5), (1.8) we have 

T,(Z, Z) = tjMP, Z)]' = ri.e^ = U (1.9) 

for every instantaneous observer Z. Hence Tj.(Z, Z) is the energy density of the particle 
flow (P, 77, m) - measured by {x,Z). Then, by the uniqueness property, the stress-energy 
tensor for a particle flow is indeed specified by (1.8). 

2 Involutive Distributions (or Foliations) 

In order to generalize the notion of particle flows, we recall the following 

Definition 3 : A p-dimensional smooth distribution D on a manifold M" is an assignment, to 
each point x G M", of a p-dimensional subspace Da; of T^M^. 

Remarks: The smoothness of D can be expressed in two equivalent ways : 

(1) Every x e M" has a neighbourhood Ux C M" on which there exists a set of smooth 
(local) vector fields {Va : a — 1, . . . ,p} such that the vectors (T4)y is a basis for the subspace 
distinguished by the distribution D^^ for every y eUx C M". Thus {Va} are said to span the 
distribution, locally. 

(2) Every x e M" has a neighbourhood Ux C M" on which there exist (n—p) independent 
smooth (local) 1-forms {9'' : k ^ p+1, . . . ,n} such that ^''Id^ = for all y eUx C M". Thus 
{9^} are called constraint 1-forms for D. If D is locally spanned by {Vq}, then 9^(ya) — 0. 
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Now, an immersed submanifold in M" is said to be an integral manifold (also called a 
'leaf') of D, if at each point x e 5^, its tangent space T^S^ coincides with the subspace of 
TxM'^. Then a distribution D is called integrable (or involutive) if through each point of M" 
there is an integral manifold of D, and the necessary and sufficient condition (Probenius's 
integrabihty condition) for D to be integrable is given by^'^ 

[K, V^] = flpV, (2.1) 

for some local functions /J^ on M". Or equivalently 

de^ = A (2.2) 

for some local 1-forms on M". An integrable distribution is called a foliation. 

3 p-Dimensional Perfect Fluids 

We recall from section 1 that a p-brane is defined by a timehke, connected, p-dimensional 
manifold immersed in a spacetime and is distinguished by a strictly positive parameter, a - 
the rest mass per unit (p — 1) -dimensional spatial volume. Thus, we introduce the following 

Definition 4 : A p-dimensional perfect fluid, (0,77,(7), in a spacetime {M'^,g) consists of 
a function rj : M ^ [0,oo), called world density, and a smooth integrable p-dimensional 
distribution D on M such that each integral manifold of D is a p-brane of rest mass per unit 
spatial volume, a. If the independent local vector fields {Vi, . . . , Vp} span D on the open set 
W C M" and a; = a{Vi A . . . AVp) / /i where det[g{Va, Vp)] = -/i^ ^ 0, so that g{uj, cu) = -a^, 
then the integral of the local number density {n — p)-form 

n^^Tjuj) (3.1) 

over a spacelike {n — p)-chain C CU, defines the total number of p-branes in C. 

Remarks: 

(1) Locally, D is spanned by a set of vector fields {V^ '■ = 1.... ,p}, which forms a basis 
for each timelike tangent space of each integral manifold (p-brane) of D. Then p-branes in 
(D, r], a) are locally characterised by a decomposable p-form on M", x = V A . . . A Vp, with 

u} = a{ViA...AVp)/i^ (3.2) 

= {ViA...AVp)/ii (3.3) 
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^;(x,x) = det[^(K,y^)] = V (3.4) 

where /i is a strictly positive, real- valued, local function. The negative sign in (3.4) reflects 
the timelike causal character of each integral manifold of D. 

(2) In (3.3), the p-form 0;^^^ (hence, u and n) is independent of the choice of the vector 
fields {Va} that span D, locally. Moreover, when restricted to ap-brane, u^P^ is the induced 
volume form on the corresponding timelike integral manifold with Q{u}p, ujp) = —1. 

(3) A p-brane of rest mass per unit spatial volume, o", is a timelike immersion with local 
parametrisation : (0, 1)^* M such that for . . . , s^) = x e M", 0*(9sa) ^ Va e T^W 
and the p-form cuo = cr{vi A ■ • • A Vp) / fj,Q satisfies 

g{ujo,uJo)^-(J^ ^0 (3.5) 

where dct[g{va,Vj3)] = —^l, and a is the Newtonian analogue of inertial energy per unit 
spatial volume of the p-brane^°. We also recall that in a curved spacetime, a p-brane is 
self-gravitating if (p is an extremal immersion (and hence, the mean curvature of vanishes). 

Thus given a foliation of a spacetime (M.g) determined by a distribution D, the correspond- 
ing multidimensioal fluid, (D, rj, a), is locally characterised by the decomposable p-form uj on 
M. In order to motivate the deflnition, (3.15), of a stress-energy tensor T for such systems, 
we now give an approximate local analysis to obtain the energy density Uz with respect to 
any instantaneous observer (x, Z) for all a; e M. 

Energy Density Uz for (D, 77,(7) : 

For any x e M", let Ux be the neighborhood of x where (0,77,(7) is locally represented 
by {00,7]), and consider ap-brane through x. Then, given an instantaneous observer {x, Z) e 
T^M^, UJ admits the following orthogonal decomposition with respect to {x, Z): 

Z A {-iz<^x) + ^z<^x (3.6) 

where 11^ = 1 + Z A is the projection operator, HzUJx is supported on {x, Z)'s rest-space, 
Z^ C TxM^^ in the sense that [IlzuJx]{Z) = 0, and Q{Z A {—lzuJx),^z^x) = 0. Now, for a 
p-brane through x, its energy per unit spatial volume, Sz, with respect to {x, Z) is defined by 

Sz = [Q{i'ZUJx,i'ZCOx)]^ (3.7) 
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To motivate this definition, we look at the relevant properties of (p — l)-form © = lzuJx — 

{cr/lla:)l'ZXx- 



(a) Prom the equations (3.2) and (3.4) we have 



p 



■X 



.a 



(3.8) 



a=l 



where the (p — l)-form is defined by 



= A . . . A K-i A K+i A . . . A Vp 



(3.9) 



Note that one of the local vector fields, say Vi, must be causal since the set {Va)x 

forms a basis for a timehke subspace of T^M. Then, g{Vi,Z) ^ since Z is timehke, and 
hence by (3.8) © 7^ 0. Since, i^© = 0Jx{Z-, Z) — and is decomposable, the dimension of 
the characteristic subspace^ of the non-zero (p — l)-form © is (n — p + 1), and hence © is 
also decomposable and can be written as 



where {Yi, . . . , l^-i} are linearly independent vectors in T^M. Now, lzQ — and the hnear 
independence of the Y^s imply 



From (3.11) it follows Ya € Z-^, and hence each Ya is spacelike and ^(0, 0) > 0. 

(b) For any normal field where g{N, Va) = for all a, we find lnO = I'NI'Z^ — 0. Then 
by (3.10) it follows that giYa, N) = and hence each spacelike Ya also belongs to - the 
tangent space of a p-brane through x G M, which is spanned by the set of independent 
vectors {Vi^, . . . ,Vpx}. 

(c) From the above characterisation of © by (3.8), (3.10) and (3.11) it is now clear that the 
non-zero (p — l)-form © is constructed from a set of linearly independent (p — 1) spacelike 
vectors, {Ya G Z-^ : a=l,...,p — 1}, which also belong to the (Lorentzian) tangent space 

of a (timehke) p-brane through x e M. Thus in Z-^ C T^M^, the (p— l)-plane formed by 
{Ya} [and hence lzuJx = {o-/ fJ'x)^-zXx] represents the 'spatial extension' of a p-brane through 
X E M with respect to an observer (x, Z) and it follows that 



e = {a/ii,){Y,A...AYp_,) 



(3.10) 



g{Ya,Z)^0 Va = l,...,p-1 



(3.11) 



Az = [G{tzXx,tzXx)]'^ > 
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is the spatial volume of the p-brane in Z^. Now, projecting the foliating p-form x oiito 
[as in (3.6)] and using ^(x, x) = — //^ [(3-4)] we also have 
= -(^2)2 + (g^)2. where 

Qz = [e?(nzXx,nzx.)]^>o 

is the volume of the p-plane spanned by {V^}, when projected into Z^ . From the above 
relation connecting /x^;, Az and Qz we may define 

7(^) ^ = [1 - (gzMz)r ^ 

and from (3.10) compute 

[^;(e,e)]5 = {aiix^)Az 

= a[i-(QzMz)r^ 

If an observer Z belongs to the tangent space of a p-brane, then Hz^^x = 0, and hence 
Qz = and ^{Z) = 1. In this case we have 



[^(0,9)] 



= a 



where © = lz'jJx- Thus £z = [Q{i-z^x-, i-z^x)]^ is indeed the energy per unit spatial volume 
of a p-brane with respect to any instantaneous observer (x, Z). 

Now, given an instantaneous observer (x, Z), we can find a 'sufficiently small' neighbor- 
hood of X e M [where curvature tensor is negligible] which (by exponential map) can be 
regarded^ as part of T^M'^. In such a neighborhood of x, we compute the 'number density' 
with respect to {x,Z) by projecting [(3.1)] into A"~^(Z-'-) [as in (1.5)]: 

nz = nz[*(77a;)]^ = 77^ ★ [Z A {-lzoJx)] (3.12) 

where we used the identity Z M<^! = (-1)''"^ ^ [iz"^] for ^ G A^iT^M"-). Then, with 
respect to {x, Z), the number of p-brancs of D intercepted by unit volume of an (n— p)-plane 
in T^M^ - orthogonal to the p-plane represented by the nonzero p-form Z A {—iz^x) ~ is 
approximately given by [as in (1.6)] 

||nz|| = [^(n^, nz)]5 = 7]x[Q{izi^x, i^z^^x)]^ (3.13) 
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Finally, taking the product of Hn^H in (3.13) and Ez [(3.7)] we find the energy density IJz 
(that (x, Z) measures) of the fluid (D, r^, a) locally characterised by oj : 

Uz = VxQ{t^zi^x, I'Zi^x) > (3.14) 



Stress Tensor for (D, 7^,17) : 

Now Uz is supposed to be equal to T{Z, Z) for every observer Z (see our discussion before 
Definition 1) for any given form of the stress tensor T. Then the structure of Uz in (3.14) 
suggests the following definition of the stress tensor T for a multidimensional fluid (D, 77, a) 
with the local representation 



where {e"} are the local basis l-forms (on M) dual to {Xi,} such that e"'{Xb) = 6"'b for 
a, b — 1, . . . ,n and La = iXa- It is clear that T is symmetric, and for any observer {x, Z) 



Then by continuity Ta:(W^, W) > for all causal W G T^M and hence T is a stress tensor. 
Furthermore, the equations (3.14) and (3.16) show that the energy density Uz is equal to 
Tx{Z,Z) for every instantaneous observer {x,Z). Hence by the uniqueness property, the 
stress tensor for the fluid, (0,1], a), is specified by the equation (3.15). We also remark that 
replacing uj in (3.15) by the energy- momentum 1-form, P, reproduces the stress tensor for 
the particle flows. 

The stress tensor T defined in (3.15) for the fluid (D, 77, a) can be written in a form which 
is more suggestive as well as convenient for applications. From (3.4) we have 



T = r]Q{LaUJ, LbUj)e"- ® 



(3.15) 



VxOil^Zi^x, l^Z^x) > 



(3.16) 



dei[g{V^,Vp)] = g{V,A...AVp, ViA...AVp) 



(3.17) 



Now, expanding iaUJ in (3.15) in terms of [(3.9)], 



p 



(3.18) 



a=l 




(3.19) 
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where = {—l)°''^^G{x°':X^) is the cofactor of the matrix element = [giVa^Vp)] and 
from (3.17), the inverse of is given by g"'^ = C"^/(— /x^). Prom these definitions and 
(3.19), we have 

T = -(a^nyg'^f'Va ® (3.20) 

Defining 

g = r^Va ® V0 (3.21) 

we note that g is simply a rank-2 symmetric tensor field on the Lorentzian manifold (M", g) 
and constructed only from the foliating vector fields {Va}. Now computing the components 
of ^ on a leaf Lp (whose tangent space is spanned by {14}), we find from (3.21) 

g{Vx,V,)=gx. (3.22) 

where we used the fact [see the definitions below (3.19)] that is the inverse of the matrix 
element gap = giVa, V^). g"'^ exists since by (3.17) g is non-degenerate : 

det[^(K, Vn)] = det[^«^] = -/xV (3.23) 

Then from (3.22) and (3.23), restriction of the symmetric tensor field g [(3.21)] onto each 
leaf Lp defines a metric on Lp - induced by the Lorentzian metric g on M. Since the leaves of 
the foliation are connected and timelike, g - restricted to a leaf - is also a Lorentzian metric 
of constant index. Thus we have a simple interpretation of (3.20) that T is proportional to 
a metric g on each integral submanifolds, and from (3.20)-(3.21) 

T = -pg (3.24) 

where the positive function (on M), p = a^rj, is the energy density measured by all observers 
tangential to the leaves Lp. 

As we mentioned earlier [see remark(2) below Definition 3 in section 2], a p-foliation of 
an n-dimensional manifold may also be prescribed hj {n — p) = q constraint 1-forms {6^} 
where 9'^{Va) = and by suitable linear combinations from the linearly independent set {6'*} 
we can get an orthonormal set of q normal fields {Nk} such that 

giNi,Nj) = 5ij 
g{N„Va) = 

Then the sets {Va} and {Nk} together form a local basis for the tangent spaces of M and 
{Nk} are the new constraint 1-forms satisfying the integrability condition (2.2). In terms of 
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these normal fields (3.24) can be written as 



T ^-p{g-Y.Nk^Nk) ; q = (n-p) 



(3.25) 



k=l 



4 Dynamics and Symnnetries of Multidimensional Fluids 

The stress tensor, T [in equivalent forms (3.15), (3.24), (3.25)], associated with a p- 
dimensional fiuid is simply a symmetric tensor field on {M^,g). However, if the spacetime 
admits Killing vector fields, K, then T may acquire new symmetries through the Einstein 
field equation 



the Ricci tensor Ric and the scalar curvature R. Since Killing vector fields generate local 
isometrics of {M,g), we have^ LkRIc = = LkR where Lk is the Lie derivation with respect 
to K. Then from the definition of the Einstein tensor G it also follows LkG = 0, and hence 



Proposition 1: If K is a Killing vector field and (0,1], a) is a p-dimensional fluid on [M^,g), 
then the world density function rj satisfles Lk?] = 0. 

Proof : Taking trace of both sides of (4.1) with T given by (3.25), we find 



where n and p are the dimensions of M and the distribution D, respectively. Since LkR = 0, 
and p = a'^Tj, it follows from (4.2) that Lk?7 = 0. □ 

Corollary 1: If T in (4.1) is the fiuid stress tensor T = —pg [(3.24)], and K is a KiUing vector 
field on (M, g), then Lk^ = 0. 

Proof : Since LkT = 0, it follows that — (Lkp)^ — p(Lk^) = 0. Then the corollary follows from 
the fact that LkP = 0, by Proposition 1, and p 7^ 0. □ 

We now discuss the significance of Corollary 1 which suggests that the symmetries of a 
spacetime (M", g) are also the symmetries of the leaves {Lp\ of a given fohation of M. First, 
we prove a consequence of Lk^ = 0, where g — g — '}l\=i ® ^k [(3-25)], the set {N^} is 
normal to the foliating vector fields {K;}a=i, gi^h ^j) — ^ij and q= {n — p). 



where G = Ric 




by (4.1), LkT = 0. 



I n 

^=y^2 



1)R 



(4.2) 
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Lemma 1 : If Lk^ = where K is Killing, then LkA^^ = ^ a^N}. where the a'^'s are real 
numbers. 



Proof : Taking the Lie derivative of g we have \-K{g — Sf=i -^i ® Ni) — 0. Since K is KiUing 
(and hence Lk^^ — 0), 

^(UiV, (^Ni + Ni® UiV,) = (4.3) 

i=l 

Now, evaluating the symmetric tensor in (4.3) on {Va, Nj} 

= ^(?(Lk7V,,K)5,,- = ^(LkA/,-,K) Vj> 

1=1 

Similarly, evaluating (4.3) on {Nj, Ni} we find 

1=1 1=1 
= g{UNuN^)+g{UN^,Ni) 

Now, substituting i — j — we have ^'(Lk, Nj, N^) Thus, LkA^j is normal to Nj as well 
as Va "ioi, and hence, expanding LkA^^ in the basis {Vi, ■ ■ ■ ,Vp, Ni, . . . , Nq\, it follows that 

Proposition 2 : If K is a Killing vector field on {M,g), then L^ou — 0, where u = a{Vi A 
. . . AVp)/n - defined by the fohating vector fields {Va} - is the local representation of the 
distribution D characterising a p-dimensional fluid {D,r],a) in a spacetime {M'^,g), and 
-fj,'' = det[g{Va,V^)]. 

Proof: In terms of the spacehke orthonormal set {Nk} 

u ^ ai-iy^P"^ {Ni A . . . A Ng) (4.4) 

where -k is the Hodge operator induced by the spacetime metric g and g{Nk, Va) = V/c = 
1, . . . , g ; Va = 1, . . . ,p. Since for any Killing field K, Lk commutes with ★ and metric dual 
operation 

Lk*(7Vi A...A7V,) = ^[LK(iVi A... A7V,)] 
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where 11^ = TVi A . . . A Nj_i A 7V,+i ... A Ng. Now, using LkNj = ^ a'^Nk from Lemma 1, it 

follows that each term in the above expansion vanishes. Hence Lk^ {Ni A ... A Ng) = 0, and 
by Lie derivation of (4.4) with respect to K we conclude that Lkcc = 0. □ 

Now we derive the dynamical consequences of (4.1) with the foliating stress tensor [(3.25)]. 
First, we recall that the Einstein tensor G is divergence-free, V ■ G = 0, and hence for any 
stress tensor (4.1) implies V ■ T = 0. In our case of interest, T is the stress tensor for a 
p-foliation, given by [(3.25)] T = —p{g — J2k=i ® -^fc)- Then, we have 

V • t = -dp + X:{V • {pNk)}Nk + X: pVN.Nk = (4.5) 

k=l k=l 

Proposition 3: For each integral manifold (p-brane) of the distribution, D, characterising a 
p-dimensional fluid (D, rj, a) in (M"^, g), the mean curvature field H = 0. 

Proof : The proposition involves a local assertion. Since the distribution, D, is integrable, 
for every point x G M"^ there exists an integral manifold, S^, passing through x, and (by 
Frobenius Theorem) there is an open neighborhood of x, U G M", where we may choose 
an orthonormal moving frame Xi, . . . , Xp, Ni, . . . , Ng such that {Xa} are tangent to 5*' and 
{Nj} are normal to S^. Since is timelike, g{Xa, Xa) — = il- Then the mean cuvature^^ 
vector field H of C has the following local representation: 

H = i:i:ea^(Vx„X„, N,)Nj (4.6) 

j=l a=l 

From (4.5), gf(V • T, Nj) = 0, which implies 

-iV, (p) + V • (pN,) +pj2 5(VA.,iVfc, N,) = (4.7) 

k=l 

Now, inserting the following expansion 

V-{pN,) = Nj{p) + pV-N, 

= Nj{p) + p X: eag(Vx.Nj, Xa)+pf: g(V^,Nj, N,) 

a=l k=l 

= Nj{p) -pfl eag{Vx.Xa, Nj) - p gi^N.Nk, Nj) 

a=l k=l 

in the equation (4.7) we have J2a=i^ag{^XaXa: ^j) — 0. Then, by (4.6), it follows that 
g{H, Nj) = for j = 1, . . . , Hence H = 0. □ 



15 



Remark: By the above proposition, the equation of motion of a p-brane, <S, in a multidimen- 
sional fluid (D,77,(t) is given by H = 0. Also, it can be shown^^ that H and oj (the local 
representation of D) are related by i^dijj\s — — 5'(H, N)\s-, for every vector field, A^, normal 
to the p-brane. Then H = implies iNdio\s — 0. 

In order to derive a further dynamical consequence of (4.5), we need the following 

Lemma 2 : The world-density function, r/, satisfies lyidri -|- t^A) = 0, where V is any vector 
field tangent to the integral submanifolds of the fiuid, (D, 77, cr) and A is some 1-form. 

Proof : From (4.5), gf(V • T, l^) = 0, which implies 

y{p) = pi:^(Viv,A^.,T/) 
fe=i 

= pj2dNk{Nk,V) 

k=l 
k=l 

where we have used the identity VAr^.A^'^ = L^^dNk and the fact that the constraint 1- 
forms {Nk} describing the p-foliation must satisfy the integrability conditions [(2.2)] dNk — 

Sj=i A-^jfc A Nj, Aj. being suitable 1-forms. Since p = a^rj, defining A = ^ A^, we have 

fe=i 

[dri + riX\{V) = 0. □ 

Proposition 4 : The world density function 77 on (M, g) satisfies d ★ (rjou) — 0. 

Proof : Prom the expression for the foliation p-form u [(4.4)] 

di.{r)uj) = a{-lY+P''{r)d{NiA...ANq) 
+dr] A {Ni A . . . A Ng)} 

Using the integrability conditions (2.2) for {Nk} we compute 

q 

d{NiA...ANg) = ^(-l)^^^d7Vfc A$*^ 

k=l 

= (E'^Ma(^iA...A7V,) 

fe=i 

where = TVi A . . . A Nk-i A Nk+i ■ ■ ■ A Nq. Then we have 

d ^ {Tjiv) = ai-iy-^P'^idT] + 7]Y. ^'k) A{NiA...ANg) 

k=l 
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Since the 1-form (dr) + r)J2 ^''k) does not have any tangential components (by lemma 2) and 
its {7Vjfc}-components do not contribute in the above equation, (rjuj) = 0. □ 

Remark : Using the coderivative operator S — (— -k d-k on differential /c-forms in a 
Lorentzian manifold (M", g), we can write d -k (rju;) = as S{r]u;) — 0. 

Thus our program - of investigating a p-dimensional fluids as a source of spacetime 
curvature - would be to solve the Einstien equations (4.1) with the stress tensor T in (3.15), 
and to specify g and {uj,rj) - the local representation of the fluid. In the next section we 
shall work out a complete solution for a fluid characterised by a 2-dimensional distribution 
on a spacetime {M^,g). 

5 Spherically Symmetric 2-Foliation 

As an application of our results in the previous sections we consider a 2-foliation (due to 
string world-sheets) of a static spherically symmetric spacetime {M'^,g) where, in the local 
chart {t, r, 9, 4>), metric g is of the form 

g = -h?{r)dt ®dt + f(r)dr ® dr + r^{de ® d9 + sin^ Odcp (g) #) (5.1) 

Using an orthonormal basis (5.1) can be written as 

g, = -e° (g) e° + + + (5.2) 

where coframes are 

e^^h{r)dt] ^ f(r)dr; ^ rdO; e^ = rsine# (5.3) 
with the dual basis given by 

Xo = {l/h)dt; X, = {l/f)dr; X2 = {l/r)de; = {1/r sin 6)8^ (5.4) 
such that e^{Xi,) — 5^b. It is clear from (5.1) that (M^,g') has four KiUing vector fields : 

Ko = dt 

Ki = sin (f)dg + cot 9 cos (f)d^ 
K2 — — cos (j)d0 + cot 9 sin 

K3 = (5.5) 
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Foliation 2-Fornn : If the spacetime specified by (5.1) is to be foliated by string world-sheets, 
we must find the appropriate 2- form oj e A^(M) satisfying [Proposition 2] 

Wjjj = (5.6) 

where the KiUing vector fields {Ki} are given in (5.5). The most general 2-form satisfying 
(5.6), on the chosen spacetime [(5.1)], must be 

u = cx{T)dt Ndr^ C2(r) sin 61^61 A d0 (5.7) 

Locally, uj is required to be decomposable, and for timclike foliation uj must satisfy Q[yj^ uS) < 
0. These two constraints together with (5.7) uniquely (up to a scalar function) specify the 
structure of the foliation 2-form so that 

UJ ^ ci{r)dt A dr (5.8) 

Without any loss of generality we can normahse (5.8) by G{u;,u;) — —1, and the foliation 
2-form is then given by 

u = h(r)f(r)dt A = e° A (5.9) 

where we used the equations (5.1) — (5.4). It is now easy to see that the contraint form ★a; 
satisfies the integr ability condition [(2.3)]: 

= -d(e^ Ae^) = (2/r)(irA*a; (5.10) 

Thus UJ, indeed, determines a 2-foliation. Now, introducing string rest-mass per unit length, 
cr, the foliation 2-form for the extended particle flow {u!,r]) is written as 

a; = (7e°Ae^ (5.11) 



Stress Tensor : From (3.15) and (5.11), the associated stress tensor is given by 

T = 7]Q{ia^,Li,uj)e"' ® 

= r]a^{g{e\ e^)e° ® e° + g{-e\ -e°)e^ ® e^} 

= 77(j2{e° e° - e^} (5.12) 

where T is expanded in the orthonormal basis given in (5.3) and (5.4). The density function, 
r], can be obtained from Proposition 4 and (5.11): 

= d ★ {r]uj) = -ad{rir^) A sin 9d9 A dcf) (5.13) 
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Then (5.13) implies dr{r]r'^) — and hence rjr'^ is constant. Thus the density function for 
2-fohation is given by 

T] = c/r^ (5.14) 
where c is some positive constant since r] is defined to be positive. 

Solution to Einstein's Equations : For complete specification of the string-field flow we must 
find the functions h{r) and /(r) from the Einstein equation (4.1) with the stress tensor in 
(5.12). For convenience (4.1) is written in the following form : 

Pa = Tabe' - {r/2)gabe' (5.15) 

where Pa = R\c{Xa, Xb)e^ are the Ricci 1-forms, and F = traceT = —2ria'^ by (5.12). Now, 
computing the Ricci forms with respect to an orthonormal basis[(5.3), (5.4)], we find 

Po = {l/fW/h)-{h'/h){f/f) + {2/rW/h)]e' 

P, = -il/f)[{h"/h)-{h'/h){f'/f)-{2/r){f'/f)]e' 

P2 = [{yrf){-{h'/h) + {f/f)} + {iy){l-{l/f)}]e' 

P3 = [{yrf){-{h'/h) + {f/f)} + {l/r'){l-{l/me' (5.16) 

From (5.2), (5.12), (5.15) and defining p = rja^, we also have 

Po = 0; Pi = 0; P2 = pe2; P3 = pe^ (5.17) 
Then (5.16) and (5.17) imply 

{h"/h) - {h'/h){f'/f) + {2/r){h'/h) = (5.18) 

{h"/h) - {h'/h){f'/f) - (2/r)(/7/) = (5.19) 

{l/rn{-{h'lh) + (/'//)} + (l/r^){l - (1/f )} = P (5.20) 
To solve these equations, first, we note that subtracting (5.19) from (5.18) gives 

{h'/h) + (/'//) = (5.21) 
Integrating (5.21) and choosing the integration constant to be 0, we have 

hf = 1 (5.22) 
Now inserting (5.21) in (5.18) and defining a(r) = /i^(r) we get 

a" + (2/r)a' = (5.23) 
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The general solution to (5.23) is found to be 

a(r) = h'^{r) - (/9 - 2m/r) (5.24) 

where f3 and m are constants with m > 0. Then using (5.21), (5.22) and (5.24) in the 
equation (5.20) we find 

(l/0(l-/5)=p (5.25) 

It is clear from (5.25) that 

(3^1^ p^O (5.26) 

Hence for non- vanishing string field fiow (cu, r/), the constant /5 can not be equal to 1. Com- 
paring (5.25) with (5.14) and using p = rja"^ we also have 

/3 = 1 - c(7^ (5.27) 

Collecting our results in (5.22), (5.24) and (5.27) the spacetime metric [(5.1)] is now given 

by 

g — —{l — ca^ — 2m/r)dt®dt-\-{l — ca^ — 2m/r)~^dr®dr 

W{de ®de + sin^ edcp ® defy (5.28) 

and from (5.22) the foliating 2-form a;[(5.11)] takes the following form : 

u^adtA dr (5.29) 

Thus the equations (5.28), (5.29) and (5.14) completely determine the local flow of a 2- 
dimensional fluid generated by radial strings in a static spherically symmetric spacetime. The 
metric in (5.28) may be interpreted as representing a spacetime associated with a particle of 
mass m (at r = 0) surrounded by spherically symmetric distribution of strings with density 
rj^c/r'^ [(5.14)]. 

Properties of the Solution (5.28) : 

(a) For ccr^ < 1, this solution has a horizon of radius 

ro = 2m/(l - ca^) (5.30) 

The equation (5.30) shows that the Schwarzschild radius for the mass m is enhanced by a 
factor (1 - c(7^)~^ > 1. 
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(b) For m — 0, there is no horizon but the spacetime has a naked singularity at r = 0. 
To see this, we substitute p — 2 (dimension of fohation) and n — A (spacetime dimension) 
in (4.2) to find the Ricci (scalar) curvature : 

R = 2ca^/r^ (5.31) 

We also compute an invariant scalar constructed from the curvature 2-forms Ra^, where 
a, 6 = 0, 1, 2, 3 : 

★(R'^^ A ^Rab) = 24mVr'^ + 8m{ca^) jr" + 2(cV^) /r^ (5.32) 

Besides the existence of the singularity at r = 0, (5.31) and (5.32) also imply that the 
spacetime remains curved with m = - that is - with the foliating strings alone, and the 
metric, (5.28), remains well-behaved even with m = and ca^ — 1. 

(c) Using (5.2), if we rewrite the foliation stress tensor, (5.12), in the form, T = —p{g — 
® — e" ® e'\ then for any causal vector V we have 

l{y,V) = -p{g{V,V)-[g{X2,V)f-[g{X,,V)f} 

> -pg{V, V) = {l/2){traceT)g{V, V) (5.33) 

where we have used the inequalities p > and g{V,V) < 0. The equation (5.33) shows that 
the stress tensor (5.12) satisfies the strong energy condition, and this also means that the 
gravitational field - generated by the fohating string world-sheets - is attractive. 

Thus the fiuid of string world-sheets gives rise to a non-trivial solution (to Einstein's 
equations) - which is non-fiat, static and spherically symmetric with a naked singularity. 
Such a network of line-like objects could be used to model a multilayer star where the 
constituents of each layer follows different equation of state. 

Modelling Star with 2-Foliation : 

Following a suggestion in reference [12], we consider a two-layer star in which the core 
consists of a (spatially isotropic) perfect fiuid, and the exterior is formed by a spherically 
symmetric distribution of strings as in (5.29). Thus the stress tensor for the core is given by 

Tc = (Pc + ®u + vgc (5.34) 

where u is a unit timelike vector field - called fiow vector field, pc is the density of the 
core and taken to be a constant, and v is the spatially isotropic pressure function. Then, 
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in the chart {t,r,9,(l)), the metric tensor Qc for the core of radius Tc is given by a special 
case (pc = constant) of the well-known Oppenheimer-Volkov^^ solution to the Einstein field 
equations : 

gc = -{l/4)[3{l-^Pcrl)-^-{l-^p,r^y^]^dt(^dt 

+ (1 - ^PcT^y^dr <»dr + r'^{d0 <^ d0 + sin^ 0d(l) ® dcj)) (5.35) 

where r e (0, Tc) and w = 9^ in (5.34). Futhermore, the isotropic pressure i/(r) can be 
obtained from the Oppenheimer-Volkov^^ equation : 

{3iy + p,f/{u + p,f = [(31/0 + Pcf/i.^0 + Pcf]{l - \pcr'') (5.36) 
where vq is the pressure at r = 0. 

The metric for the spacetime region foliated by the radial strings [(5.29)] is taken as 
[(5.28)] : 

Qs — —{1 — q — 2m/r)dt <S> dt + {1 — q — 2m/r)~^dr <S> dr 

+r^{d9 ®de + sin^ Odcl) ® d(j)) (5.37) 

where q = ca^ [from (5.28)] is a positive constant. In (5.37) we require r G (r^r^) with 
as star-radius, and > Tg. For r > r^, we have the Schwarzschild vacuum metric : 

= —{l — 2m-v/r)dt<Sidt+{l — 2mv/r)~^dr<Sidr 

+r^{d9 ®de + sin^ ed(p ® d(p) (5.38) 

To complete our model we need to match these metrics continuously across the boundaries 
of different layers. The matching condition for Qs and Qc r — Tc is given by Qs 
and hence from (5.35) and (5.37) we have 



= 9c 



2m = r,{^-p,Tl - q) (5.39) 
Similarly, from (5.37) and (5.38), the continuity of gs and gi^ at r = implies 

2m^ = 2m -f qVs = ^Pc^l + q{rs - Tc) (5.40) 
where > Tc. Now, if we impose the condition [see (5.35)] 

IpctI < 1 (5.41) 
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then it follows from (5.39), (5.40) and (5.41) that g < 1, and 



Tc > 2m/(l -q) ; > 2m^ (5-42) 

The above conditions, (5.41)-(5.42), ensure that the spacetime regions specified by the cor- 
responding metrics [(5.35), (5.37) and (5.38)] are static and free from any singularities. We 
also remark that (5.41) naturally follows from the condition that prohibits gravitational col- 
lapse of the core. To see this, first we note that the pressure = at the core-surface r = Tc, 
and this implies [from (5.36)] 

1 = [(3^.0 + Pc)7(^o + Pom - \pcrl) (5.43) 
From the above equation (5.43), we find 

^Pcrl = iuo{2uo + Pc)/(3i/o + Pcf (5.44) 

The right side of (5.44) can be easily seen to be an increasing function of the central pressure 
uq. However, evaluating the limit of (5.44) as i^o we get 
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^Pcrc 



i^O— >c» 9 



(5.45) 



The equation (5.45) shows that there exists a maximum 
Tc with the given density Pc, and hence 



for the core-radius 

i/Q— >00 



IpctI < \pcRl < 1 (5.46) 

The above inequality provides the validity of the condition (5.41) which leads to the model 
of a non- collapsing star with spherically symmetric distribution of strings. 

Our discussion on 2-foliation suggests that a complete general-relativistic theory of mul- 
tidimensional perfect fluids may enable us to investigate a new class of collapse problems, 
the possible formation of horizons and the nature of the associated singularities. 

6 Conclusion : 



We emphasize that the results in this paper strictly follow from the concepts intro- 
duced in Definition 4 which offer a precise description of p-dimensional fluids. Consequently, 
the foliation p-form cu [with Q{u!,u!) < 0] together with the uniqueness property of stress- 
energy tensors directly lead to our formulation of the dynamics and symmetries of such 
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self-gravitating systems. Then the local decomposable p-form to which defines a timelike 
p-foliation of a spacetime and the world density function 77, enable us to introduce the local 
number-density [(3.1)] of the leaves of the fohation. They also give rise to the concept of 
spatial volume (and energy) of ap-brane with respect to any observer [remarks (a), (b) and 
(c) after (3.7)]. These two ingredients then naturally motivates a precise definition of the 
stress-energy tensor T[(3.15)] for multidimensional perfect fiuids. 

It is interesting to observe, from the equation [(4.2)] relating trace(T) and the the scalar 
curvature R, that a 2-dimensional spacetime (n = 2) does not admit any massive particle 
(p = 1) fiow, however, only massless fiows are consistent with the field equations. In fact, a 
slight modification of our Definition 4 permits construction of null foliations, where u; satisfies 
Q{u;,u;) — 0. 

Furthermore, u; provides the local description of the dynamics [Proposition 3 and Propo- 
sition 4] of p-branes and carries the symmetries [Proposition 1 and Proposition 2] of a fohated 
spacetime, and we have demonstrated, in section 5, that the spacetime symmetries determine 
to a large extent the structure of fohations [(5.7)] and hence the associated stress tensor. 
This example suggests the posssibihty of other non-trivial solutions to the Einstein equations 
in a spacetime region dominated by extended objects like strings and membranes. 

Appendix : Definition of Q 

Given a semi-Ricmannian manifold {M"',g), the symmetric bilinear form Q on the vector 
space AP(M) of differential p-forms on M is, first, defined^ on decomposable p-forms and then 
the definition is extended to any p-forms by linearity. Consider two decomposable p-forms 
6 and $ given by 

© = cci A . . . A CKp ; $ = /?i A . . . A /9p 

where ctj. Pi e A^{M). Then, by definition, Q{Q, $) = det[g{ai, Pj)] - where a is the metric 
dual of the 1-form a. Now, for any pair of p-forms cu and x, Q has the following useful 
properties : 

a;A*x = e?(^,x)*i ; Gi^,*x) ^ i-iTGi^^,x) 

where -kl is the volume form on M induced by the metric g, and s is the index of g. 
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